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Abstract 

An Improved method Is presented for cou- 
pling a boundary layer code with an unsteady 
Invlscld transonic computer code In a quasi - 
steady fashion. At each fixed time step, the 
boundary layer and invlscld equations are suc- 
cessively solved until the process converges. 
An explicit coupling of the equations Is des- 
cribed which greatly accelerates the convergence 
process. Computer times for converged viscous- 
in viscid solutions are about 1.8 times the com- 
parable invlscld values. Comparison of the 

results obtained with experimental data on three 
airfoils are presented. These comparisons 

demonstrate that the explicitly coupled 
vlscous-lnviscid solutions can provide efficient 
predictions of pressure distributions and lift 
for unsteady two-dimensional transonic flows. 

Nomencl ature 

CT computational test case 

C z first harmonic lift coefficient 

a due to pitch 

C ma first harmonic pitching moment 

coefficient due to pitch 
Cp pressure coefficient 

c|: entrainment coefficient 

c airfoil chord 

F airfoil surface function 

F ! ,F 2 boundary layer function, eq. (5) 

TT boundary layer shape factor 

k reduced frequency, wc/2U 

M free stream Mach number 

t time 

U free stream velocity 

x streamwise coordinate relative to 

leading edge 

y coordinate normal to freestream, 

positive up 

A( ...) Indicates junp In ... 

oq mean angle of attack 

ai dynamic pitch angle 

Y ratio of specific heats 

Y* Y* = 2 - (2 - y)M 2 

6 airfoil thickness ratio 

6* boundary layer displacement thickness 

6 boundary layer momentum thickness 

4 > Invlscid perturbation velocity 

potential 

4 >V perturbation velocity potential with 

viscous interaction 
u) angular frequency 

Introduction 

The Importance of viscous effects in accu- 
rate predictions of steady and unsteady tran- 
sonic aerodynamic loading Is well-known. 1 While 
Navier-Stokes computer codes can provide solu- 
tions 2 which Include viscosity, their cost Is 


prohibitive for routine use. As a result, 
extensive efforts are underway by various 
researchers to account for unsteady viscous 
effects by coupling a viscous boundary layer 
model with an otherwise invlscld analysis. 6 
As commonly implemented, the Invlscid outer flow 
solution provides the surface pressure distribu- 
tion which is needed to solve the boundary layer 
equations. This yields the boundary layer 
thickness distribution which is used to modify 
the airfoil surface tangency boundary condition 
for the next iteration of the outer Invlscld 
flow solution. 

For steady flow problems. It has been 
demonstrated that this direct solution technique 
converges to consistent solutions of the boun- 
dary layer and inviscld equations which show 
good agreement with experiments 7 . These con- 
verged interactive solutions are referred to as 
self-consistent solutions. For unsteady flow 
problems, these iterative solution techniques 
are impractical due to the hlcji computational 
expense associated with a large number of Itera- 
tions. 

Rizzetta 1 coupled the steady Integral lag 
entrainment boundary layer model of Green 8 with 
the LTRAN2 9 unsteady transonic code in a quasi - 
steady manner using a non-iterative implicit 
technique. Guruswamy 3 applied the method of 
Rizzetta in a study of viscous effects on oscil- 
lating airfoils with the result that up to 8000 
time steps per cycle were required to obtain 
reasonably accurate answers. For step sizes 
even smaller than this the computer code became 
unstable. An alternative procedure was reported 
by Houwink 4 who used an explicit method of vis- 
cous-i nvi scld coupling and obtained satisfac- 
torily converged solutions with 120 time steps 
per cycle and computer times only 30% higher 
than the invlscid cases. 

The present study began by incorporating 
Rlzzetta's 1 boundary layer algorithm into the 
XTRAN2L code. 10 Several changes have been made 
to this procedure, the most important being the 
ability to iterate the vlscous-lnviscid solu- 
tions at each time step and the inclusion of an 
explicit coupling procedure. The boundary layer 
equations are numerically integrated from a 
specified transition point to the downstream 
boundary. The paper discusses the modifications 
which have been made to the computational algo- 
rithm and presents comparisons of the results 
with experimental data on three airfoils. 

Analysis 

The inviscid code used in this study Is the 
XTRAN2L computer code described by Whitlow. 10 
This code is a modified version of the LTRAN2 
code developed by Ballhaus and Goorjian. 9 The 



XTRAN2L code solves the complete 2-0 
small disturbance (TSD) equation given 
tlon (1) and Includes nonreflecting 
boundary conditions. 


transonic 
by equa- 
far- field 


4k 2 M 2 4kM 2 

6 2/3 U + & Z/3 * xt 

= C(l-M 2 )/6 2/3 - M 2 ( y* +l)* x ]* xx ♦ * yy 


The disturbance velocity potential * is 

normalized by cU a 2 / 3 , c Is the airfoil chord, 
6 is the airfoil thickness ratio, and U Is the 
freestream speed. The spatial coordinates, x 

a /»l/^' anrf the i t1n ' e > t, are normalized hy c, 

c/s , and u" 1 respectively, where o> Is the 

frequency of unsteady motion. The reduced 
frequency (based on semi -chord) k = uc/?U, M 

, the ,, fre 1u2 ream u Mach number, ‘and 
Y " *( 2 " y) m » where y is the ratio of 

specific heats. 

The boundary conditions on the airfoil and 

wake for the Invlscld code are: 


airfoil : a" = F - + F~ 

y x t 


( 2 ) 


wake: A$ y = 0 

4( *x + 


(3) 


♦ t ) = 0 (4) 

where the ± refers to the airfoil upper or lower 
surface, the function F(x,t) denotes the airfoil 
surface, and A(...) indicates a jura in the 

bracketed quantity. 


To account for the viscous boundary layer in a 
2i* si ' st ^] manner, equations {?.) and (3) are modi- 
fied as follows (see Houwink 5 and Rizzetta 1 ): 


airfoil: <p~ 


F v + + (&>j 


wake: 


t 

/ 6 * , 


k 6^x 


A( *y> = a( k'x 


(Za) 

(3a) 


where 6* is the boundary layer displacement thick- 
ness determined from Green's lag entrainment equa- 
tlons. ^ 

Equation (2a) is a direct extension of the air- 
foil boundary condition + as given hv 
Rizzetta to include the term Ff which accounts 
for the time dependence of the airfoil motion in the 
boundary conditions. Note that althou^i equation 

is a commonly used form to describe analyti- 
cally the inclusion of the boundary layer equations, 
' pr f' s ' ™" ner which this equation Is Incor- 
porated into the computational algorithm can be a 
ke »w ^ 1 ea / ture °* the efficiency and accuracy of the 
method (see ref. 11, for example). This feature Is 
discussed further in the following section. 

a f„^?/ iSp1 , aCe .T nt u th1ckness * Is computed as 
a function of the boundary layer shape factor 

FT and the momentun thickness e, which are deter- 
mined, together with the entrainment coefficient 
C E» from Green's lag entrainment equations: 


°x " f l + f 2 *xx 

ff x - f3 + U *xx 
( c £)x = f 5 + hx 


of 9 H C E> and other parameters which are 

escribed in detail in references 1 and 8, 


Computational Procedure 


.... a.yontnms for the viscous calc... 

! at ■ tionaTn* m -] ar t0 th0Se of Rl^etta', w ere 
1 ^ n eta ; s he found. Several modi! 
neat ions to the procedures of reference l 

Jode and" i 1 . ncorpor ^ed into the present compute? 
code and these modifications are enumerated In 
the following paragraphs. enumerated in 


1. In the original algorithm, the viscous 

equations are integrated using a Runge-Kutta 
algorithm given the value of (C D ) X at each 
outer flow grid point on the airfoil and wake. 
The integration spatial step size was found to 
be too coarse and a finer step size is used for 
integration of the boundary layer equations. 
Ten boundary layer grid points are used between 
each pair of outer flow grid points. The 
required values of (C p ) x at the refined 

viscous grid points are determined by 
interpolation. 

2. The empirical viscous wedge was elimi- 
nated, and integration of the boundary layer 
equations begins at a specified transition 
point. This allows a consistent description 
of the boundary layer from the specified 
transition point to the downstream boundary. 
Green s lag entrainment equations are simply 
integrated through the shock with no adverse 
effects. For results in this paper, transition 
is fixed at 10% chord unless stated otherwise. 

3. Boundary layer smoothing was introduced 

to reduce instabilities. The values of * x , 
<1>xx us ed in the boundary layer calculations as 
well as the values of 5* and 5 X * are smoothed 
by averaging the values at 3-5 adjacent qrid 
points. 3 


4. An option was added to allow for 
iterating the vi scous-i nviscid solutions at each 
time step. In Rizzetta's analysis, 1 the last 
term on the right of equation (2a) is 

Implemented as: 


(-) 


- F i + f 2 


F, + — 
1 x. 


T xx 

2F 0 


i +1 


- x 


i-1 

n+1 


L( ~^ " ^ J ) 

x l+i ‘ x i 


(5) 


<fc n+1 

,*i.J 


, n+ l 

* 1-1 


where F* and F2 are functions described in 
reference 1. Now and <^ + j In this equa- 

tion implicitly couple’the two tine levels n and 


2 


n+1 at the two streamwise stations x^ + -j and 
x-j • In the present analysis, these Implicit 
equations are Iterated at each time step. 


Figure la presents a flow chart which des- 
cribes the Iteration procedure. For t * t n +i 
the boundary layer functions F lt F 2 are calcu- 
lated using the potential solution <j> n from 
the previous time step. During the y-sweep 
those , functions are used In eq. (5) to 

Implicitly couple the boundary layer with the 
Invlscid outer flow. The resulting solution for 
the coupled potential Is denoted by <{>V n+1 . 
This coupled solution, 4>V n+1 , Is then used to 
re-calculate the boundary layer parameters and 
the process Is repeated. Most applications 
investigated have required 5-25 iterations for 
convergence. 

5. The coupling between the viscous boun- 
dary layer and the invlscid code was made expli- 
cit. This alternate approach to viscous-invis- 
cld coupling Is based upon the work of 
Houwink. 4 The explicit coupling between the 
boundary layer and the Invlscid solution on the 
airfoil (but not on the wake) Is implemented 
by the direct use of equation (2a) in the 
expression for the airfoil downwash. That Is, 
at time level t * t n +i, the last term on the 
right of equation (2a) is evaluated by using 
values at the previous time step t s t n as 
follows: 



a) ORIGINAL VISCOUS COUPLING 



b) EXPLICIT VISCOUS COUPLING 


Fig, 1 Flow charts for boundary layer coupling 
Iterations. 
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( 6 ) 


This procedure does Involve a lag of boundary 
layer displacement thickness by one time step, 
but it does not Involve an implicit coupling 
between different time levels at different 
streamwise stations. 


A flow chart which Indicates the explicit 
coupling algorithm is shown in figure lb. As 
the figure indicates, the modifications for the 
explicit viscous coupling are included during 
the y-sweep of the alternating-direction-impli- 
cit (ADI) solution. The x-sweep calculations 
are Identical for both forms of coupling. For 
computational efficiency, the self-consistent 
iterations are only done over the y-sweep. In 
practice, converged self-consistent solutions 
have been obtained in one Iteration for most 
cases, resulting in more accurate viscous-invls- 
cid solutions with reasonable computer times. 

Results and Discussion 
Convergence Studies for Original Coupling 

In order to test the validity of the self- 
consistent iteration procedure as described by 
modifications 1 through 4, extensive calcula- 
tions have been performed on several airfoils. 
The results presented In this section are for 
the NACA 64A010 airfoil (theoretical section). 
This airfoil, at a Mach number of 0.78 and one 
degree angle of attack has a moderately strong 
shock near midchord that is typical of the cases 
of interest. The results obtained for this 
airfoil are similar to the results obtained for 
the other airfoils studied. 

The calculation procedure was as follows: 
Steady results for an invlscid analysis were 
calculated and used as a starting solution for a 
steady viscous solution. The steady viscous 
solution was then used as a starting solution 
for the unsteady viscous calculations. Unless 
otherwise stated, the unsteady results presented 
herein were run at 360 time steps per cycle and 
transient effects were sufficiently damped out 
after two cycles to obtain the harmonic compon- 
ents by a Fourier analysis. The number of 
viscous iterations per time step was varied and 
various unsteady quantities such as pressure, 
lift, and moment were compared as the number of 
Iterations was increased. Comparisons with the 
results of Guruswamy 3 indicate that the effect 
of these iterations is closely equivalent to 
taking smaller time steps. That Is, 360 time 
steps per cycle with two Iterations per time 
step gives results similar to using 720 time 
steps per cycle. However, the Iterative 
procedure does converge, as will be shown below, 
whereas the use of successively smaller time 
steps can lead to stability problems 3 . 

Fig. 2 shows pressure distributions on the 
airfoil upper surface for steady flow at one 
degree angle of attack and 0.78 Mach number. As 
the figure indicates, a moderately strong shock 
is located near midchord with the viscous shock 
location slightly upstream of the Invlscid 
shock. 
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Fig. 2 Upper surface steady pressure for NACA 
64 AO 10 (theoretical) airfoil: original viscous 

coupling, oq * 1°, M ** 0.78 


The first harmonic components of unsteady 
pressure distributions for one degree oscilla- 
tion in pitch (aj* 1°) with k * 0.1 for varying 
numbers of iterations are shown in Fig. 3. For 
one iteration (Fig. 3a), the viscous solution is 
quite different from the inviscid solution. For 
five iterations (Fig. 3b), the viscous and 
inviscid solutions are nearly coincident except 
in the vicinity of the shock. Similar to 
results for the steady solution, the unsteady 
viscous solution has a somewhat weaker shock 
wave at a slightly upstream location. Compari- 
son of Figs. 3b and 3c shows that the viscous 
solution after 10 Iterations Is nearly identical 
to the solution with five iterations except for 
some minor differences in magnitude near the 
shock. This clearly demonstrates the importance 
of obtaining converged self-consistent solutions 
in assessing the effects of viscosity on pres- 
sure distributions. 

Fig. 4 shows the real and imaginary parts 
of the first harmonic component of the unsteady 
lift and moment as the number of iterations per 
time step is varied. For both lift and moment, 
the real part of the coefficient changes signi- 
ficantly as the number of iterations is 

increased from one to five, thereafter settling 
down and asymptotically approaching the con- 
verged value. The moment coefficient changes 
sign between the second and third Iterations 
before asymptotically approaching the converged 
value. Since the moment coefficient Is calcu- 
lated about the quarter-chord for this case, 
some of this variation may be due to the sensi- 
tivity of the coefficient to small changes in 
the flow. Although not shown on this figure, 

calculations have been carried out for up to 22 

iterations per time step and the results were 

identical to those obtained for 10 iterations. 
The imaginary part of the unsteady lift changes 
very little during the iterations whereas the 
trend for the imaginary part of the unsteady 

moment is quite similar to that of the real 

part. Also shown on the figure are the corre- 
sponding results for the inviscid analysis. The 
converged solution Indicates that for this case 
the effects of viscosity are very minor and pri- 
marily result In a small reduction In the magni- 
tude of both the unsteady lift and moment. 





Fig. 3 Effect of number of viscous Iterations 
for original viscous coupling method on upper 
surface unsteady pressure for NACA 64A010 (theo- 
retical) airfoil: oq * 1°, a\ a 1°, N * 0.78, 

k « 0.1. 
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b) UNSTEADY MOMENT 


Fig. 4 Effect of number of viscous Iterations 
for original viscous coupling method on unsteady 
forces for NACA 64A010 (theoretical) airfoil: 
on * 1°. ai 0.78, k « 0.1. 

Explicit Coupling Studies 

The computer costs associated with the 
large number of iterations required for 
converged self-consistent solutions with the 
original formulation of vlscous-lnviscid 
coupling is quite high. Each viscous Iteration 
Increases the computer time by about 80% of the 
inviscid solution time. A practical alternative 
for reducing the number of iterations has been 
found to be the explicit coupling between the 
inviscid equations and the boundary layer dis- 
placement thickness described previously in mod- 
ification number 5. The explicitly coupled 
equations compute a self-consistent solution 
with very few iterations, usually just one. 
Hence, explicit coupling results in converged 
self-consistent solutions in computer times 
which are about 1.8 times the comparable invis- 
cid values in many cases. 

Fig. 5 shows the first harmonic of the 
unsteady pressure distribution on the upper sur- 
face of the NACA 64A010 airfoil at one degree 
angle of attack, Mach number of 0.78, ol \ - 1°, 
and k ■ 0.1 as computed using the explicit boun- 
dary layer coupling. Comparison of Fig. 5a with 
Fig. 3c shows that the pressures obtained with 
one iteration of the explicit coupling are prac- 
tically identical to those from the original 
coupling with 10 iterations. As shown by Fig. 


5b, the results for 10 iterations with explicit 
coupling are the same as those obtained with one 
Iteration, thus demonstrating that this solution 
Is In fact a converged self-consistent solution. 

Comparison with Experiments 

In this section calculations from the pre- 
sent analysis with the explicit vlscous-lnviscid 
coupling are compared with experimental data for 
several of the computational test cases selected 
by the AGARD Structures and Materials Panel 12 as 
well as some published data for the MBB-A3 air- 
foil. 13 Additional Information on the experi- 
mental configurations may be found In the refer- 
ences. 13 

NACA 64A010 . The experimental results pre- 
sented here are for the model tested at the NASA 
Ames Research Center and are taken from Chapter 
2 of reference 12. Fig. 6 shows steady pressure 
distributions on the airfoil lower surface for 
the inviscid, viscous and experimental results. 




Fig. 5 Unsteady upper surface pressure with 
explicit viscous coupling compared with Inviscid 
results for NACA 64A010 (theoretical) airfoil: 
m 1°* ®| * 1°» k * 0.1. 
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The Mach number Is 0.796 and the angle of attack 
is -0.21 degrees. These steady results corre- 
spond to AGARD Computational Test (CT) case 
5 12 . As for the previous results, the viscous 
effects are small. The predicted shock location 
for the viscous analysis Is slightly forward of 
the Inviscid shock location and agrees well with 
the experimental result. Downstream of the 
shock, both the viscous and Inviscid analyses 
agree fairly well with the experiment. 

Fig. 7 presents unsteady pressure distribu- 
tions on the lower surface of the airfoil for 
AGARD CT case 5. The first harmonic components 
are compared with the experimental values for 
ai * 1° and k = 0.101. The viscous shock pulse 
is about two percent chord upstream of the 
Inviscid shock pulse and correlates better with 
the experimental values. 



Fig. 6 Comparison between analytical and exper- 
imental steady pressure for lower surface of 
NACA 64A010 (experimental) airfoil: Explicit 

viscous coupling, oq * -0.21°, M 3 0.796 

(ref. AGARD CT Case 5). 



Fig. 7 Comparison between analytical and exper- 
i mental unsteady pressure for lower surface of 
NACA 64A010 (experimental ) airfoil: Explicit 

viscous coupling, 1 iteration, oq * -0.21°, 
oi - l w , H - 0.796, k - 0.101 (AGARD CT Case 5). 




c) IMAGINARY PART OF UNSTEADY MOMENT 

Fig. 8 Unsteady forces versus reduced f requeue 
for NACA 64A010 (experimental) airfoil. 
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In Fig* 8 f unsteady lift and moment coef- 
ficients versus frequency are presented. As 
shown in Fig. 8a f the viscous solution for the 
lift coefficient Is closer to the experimental 
results than the Inviscld solution and, for 
reduced frequencies greater than 0.1, agrees 
well with the experimental data. The calculated 
moment coefficients shown In Figs. 8b and 8c 
have the same trends as the experimental data 
although the actual values are significantly 
different. The source of this difference is 
unknown and further investigation is needed. 
The results. do show that viscous effects are not 
the dominant effect In this difference. 

For this airfoil, over the range of reduced 
frequencies investigated, the explicit coupling 
method provides efficient viscous predictions of 
unsteady pressure distributions and lift coeffi- 
cients which agree better with experiment than 
the invlscid results. 

MBB-A3 . The experimental pressure distribu- 
tion for the supercritical MBB-A3 airfoil for 
steady flow is taken from reference 13. 

Fig. 9 presents comparisons of steady pres- 
sure distributions for the experiment, the vis- 
cous and the invlscid analyses for the design 
condition of the airfoil. Note that the calcu- 
lated values were obtained by using the actual 
experimental values for the Mach number and the 
angle of attack rather than values which have 
been adjusted to match flow conditions In the 
wind tunnel as is frequently done In comparisons 
with this particular data. For the viscous 
solution, the shock is located about four per- 
cent chord forward of the Invlscid shock loca- 
tion and is somewhat weaker. However, a large 
discrepancy between the experimental values and 
the viscous prediction still exists and further 
investigation Is required. Possible explana- 
tions for this discrepancy Include Inadequate 
theoretical modeling of the shock-boundary layer 
interaction, trailing edge effects, or inade- 
quate modeling of the flow conditions In the 
wind tunnel . 



Fig. 10 presents the unsteady pressure dis- 
tribution for the viscous and inviscld solutions 
for a reduced frequency of 0.1 and ax = 0.5°. 
The results are analogous to those for the 
steady case. On the upper surface of the air- 
foil, the shock location for the viscous solu- 
tion Is about five percent chord forward of the 
Invlscid shock location and the shock strength 
for the viscous solution is considerably weaker 
than that of the invlscid solution. On the 
lower surface of the airfoil, the unsteady 
pressure distributions for both solutions are 
essentially the same. ; 

NLR 7301. The experimental results for the 
supercritical NLR 7301 airfoil are taken from 
Chapter 4 of reference 12. Calculated results 
presented herein correspond to the test condi- 
tions for AGARD CT cases 3 and 5 12 although some 
differences do exist between these values and 
actual wind tunnel test conditions. An excel- 
lent discussion of this point is given by 
Lambourne in Chapter 0 of reference 12. 


Viscous 



Viscous 

Inviscld 



b) LOUER SURFACE PRESSURE 


Fig. 9 Comparison of analytical and experimen- 
tal steady pressure for supercritical HBB-A3 
airfoil: explicit viscous coupling, a* -1.5°, 

irr '' 


Fig. 10 Coi*>ar1son of viscous and Inviscld ana- 
lytical unsteady pressure for supercritical 
M1B-A3 airfoil: explicit viscous coupling, 

Oq «* 1.5°, aj * 0.5°, M = 0.765, k - 0.1. 
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Fig. 11 presents steady pressure distribu- 
tions for experimental and analytical cases with 
M * 0.7 and og ■ 2°. These values corre- 

spond to the AGARD CT cases 3 and 5. For the 
viscous solution, transition Is fixed at 30% 
chord to correspond to the experimental configu- 
ration 12 . The shock location for the invlscld 
solution Is about 20% chord downstream of the 
experimental shock location. The viscous solu- 
tion has a shock location and strength nearly 
coincident with the experimental result. On the 
airfoil lower surface, pressures from the vis- 
cous solution are slightly below the experimen- 
tal values although they agree better with the 
experiment than the invlscld values. This close 
agreement between the experiment and the viscous 
solution for this difficult case may be somewhat 
fortuitous due to the previously mentioned dif- 
ferences between the parameters for the AGARD 
cases and the wind tunnel tests. 




Fig. 11 Conparlson of analytical and experimen- 
tal steady pressure for supercritical NLR 7301 
airfoil: explicit viscous coupling, oq » 2°, 

N - 0.7 (ref. AGARD CT Case 3 and 5). 


In Fig. 12, unsteady pressure distributions 
are plotted corresponding to AGARD CT case 3 
with a i - 0.5° and k = 0.072. The unsteady vis- 
cous results for this airfoil were computed with 
720 time steps per cycle and two viscous itera- 
tions (explicit coupling) per time step. As 

Fig. 12a shows, on the airfoil upper surface the 
viscous solution agrees much better with the 
experiment than the invlscld solution. The vis- 
cous shock pulse is slightly downstream of the 
experimental location, significantly weaker and 
broader. As shown in figure 12b, on the airfoil 
lower surface viscous effects are small and the 
agreement between calculated results and 
experiment is good. 

Fig. 13 presents unsteady pressure distri- 
butions for a higher reduced frequency of 
k - 0.192 with cm « 0.5° (AGARD CT case 5). As 
with the previous case, the viscous solution 
agrees quite well with the experiment whereas 
the Invlscld solution Is considerably differ- 
ent. Fig. 13a shows that the location of the 



Fig. 12 Comparison of analytical and experimen- 
tal unsteady pressure for supercritical NLR 7301 
airfoil: explicit viscous coupling, 2 Itera- 

tions, ao * 2°, <n - 0.5°, H = 0.07 k * 0.072 
(AGARD CT Case 3) 




Fig. 13 Comparison of analytical and experimen- 
tal unstea<(y pressure for supercritical NLR 7301 
airfoil: explicit viscous coupling, 2 itera- 

tions, oq * 2°, ox * 0.5°, M * 0.7, k * 0.192 
(AGARD CT Case S) . 
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viscous shock pulse Is correctly predicted 
although it Is much narrower and some differ- 
ences are noted in amplitude. On the airfoil 
lower surface, viscous effects are small and the 
calculated results agree well with the experi- 
ment as shown in Fig. 13b. 

Conclusions 


This paper has presented a study of self- 
consistent solutions for viscous-inviscid inter- 
actions in unsteady two-dimensional transonic 
flow. Two different numerical techniques have 
been Investigated for coupling the quasi -steady 
viscous boundary layer with the inviscid solu- 
tion: 1) The original method of Rizzetta 1 with 
the addition of self-consistent iterations; and 
2) Explicit coupling based upon the work of 
Houwink**. Both methods are based upon the 2-D 
XTRAN2L inviscid transonic computer code coupled 
with a turbulent viscous boundary layer repre- 
sented by Green's lag entrainment equations. In 
this study, the viscous wedge in Rizzetta's 
analysis has been eliminated and the transition 
point is specified. Each method successively 
solves the boundary layer and inviscid equations 
at each fixed time step until the process con- 
verges. The results demonstrate that both 
methods converge to self-consistent solutions. 
However, method 2, which incorporates explicitly 
coupled viscous-inviscid equations, computes a 
converged self-consistent solution in one or two 
iterations whereas method 1 requires 5 to 25 
iterations to converge. The computer code with 
explicitly coupled equations yields computer 
times which are reduced by a factor of 20 from 
those of method 1. Accurate, self-consistent, 
converged, viscous-inviscid solutions are 
obtained in computer times that are about 1.8 
times the comparable inviscid values. 

Comparisons of unsteady forces and pressure 
distributions computed with the original cou- 
pling method for various numbers of iterations 
have demonstrated the importance of obtaining 
converged self-consistent solutions for an 
accurate assessment of viscous effects. 

For the NACA 64A010 (experimental ) airfoil, 
unsteady lifting forces obtained by the viscous 
calculations with explicitly coupled boundary 
conditions agree well with experimental results 
for reduced frequencies greater than 0.1 over 
the range studied. Unsteady moment coefficients 
Indicated the same trend as the experimental 
results although the magnitudes were signifi- 
cantly different. 

Comparisons of steady pressure distribu- 
tions for the supercritical MBB-A3 airfoil show 
that the viscous shock location (explicit cou- 
pling) is about 4% chord forward of the inviscid 
shock location although a large discrepancy 
still exists between the viscous results and the 
experiment. The source of this discrepancy is 
not known and further investigation is 
required. 

For the supercritical NLR 7301 airfoil, 
comparisons of experimental results with viscous 
solutions computed by the explicitly coupled 
equations clearly demonstrate that this type of 
analysis can yield accurate predictions of both 


steady and unsteady pressure distributions for 
this difficult case. Although some differences 
in shock strength do exist for the conditions 
Investigated, both steady and unsteady shock 
locations are predicted with good accuracy. 

The results presented demonstrate that 
self-consistent viscous solutions computed with 
the explicit coupling algorithm can provide 
efficient predictions of pressure distributions 
and lift for unsteady transonic flow which 
correlate better, sometimes significantly 
better, with experimental values than the invis- 
cid solutions. 
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